Spin, charge and orbital fluctuations in a multi-orbital Mott insulator 
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The two-orbital degenerate Hubbard model with distinct hopping integrals is studied by combining 
dynamical mean-field theory with quantum Monte Carlo simulations. The role of orbital fluctuations 
for the nature of the Mott transition is elucidated by examining the temperature dependence of spin, 
charge and orbital susceptibilities as well as the one-particle spectral function. We also consider the 
effect of the hybridization between the two orbitals, which is important particularly close to the Mott 
transition points. The introduction of the hybridization induces orbital fluctuations, resulting in the 
formation of a Kondo-like heavy-fermion behavior, similarly to / electron systems, but involving 
electrons in bands of comparable width. 



PACS numbers: 71.10.Fd, 71.30.+h 



I. INTRODUCTION 

Strongly correlated electron systems with multi-orbital 
bands pose a variety of intriguing problems. One 
of the recently debated topics is the orbital-selective 
Mott transition (OSMT) in highly correlated d-electron 
systems i2iiiS4& It is a fundamental issue of multi-orbital 
systems whether Mott-transitions would take place in 
sequence or simultaneously for all bands, if correlation 
would gradually be turned on. There are however also 
specific materials which have been discussed in this con- 
text such as the calcium-doped single layer strontium 
ruthenate Ca2- x Sr x Ru04£ and the ternary nickel oxide 
La n+ iNi n 03 n+ iriS where the chemical substitution (or 
the change in the temperature) may trigger the OSMT 
in the t^g (e g ) orbitals in the former (latter) case. 

The extensive studies on the Mott transition in the 
multi-orbital systems clarified that the competition be- 
tween the intra- and inter-orbital interactions as well as 
the Hund coupling plays a key role to determine the na- 
ture of the Mott transition^ It was found that under spe- 
cial conditions in a two-band system the Mott-transitions 
may merge to a single one, but would split for a generic 
form of the model. In particular, the presence of Hund 
coupling seems to be essential to observe distinct tran- 
sitions. These conclusions were drawn from the analysis 
of the quasi-particle weights computed at zero tempera- 
ture. In order to characterize the transitions the behavior 
of the spin, charge and orbital fluctuations provides ad- 
ditional valuable information. A systematic study of the 
temperature dependence of certain susceptibilities will 
give us the necessary insight to analyze in particular the 
electronic degrees of freedom which are localized through 
the Mott transition. 

The above discussions on the Mott transition are re- 
stricted so far to systems for which the bands do not 
hybridize, but are coupled to each other only through 
electron-electron interactions. However, the hybridiza- 
tion between the bands may be important in some 



compounds. 1 -- In particular, this effect could give rise to 
a qualitative change in the phase diagram, when there 
occurs the OSMT, for which the intermediate phase ap- 
pears with one orbital localized and the other itiner- 
ant. One thus naively wonders whether Kondo-like heavy 
fermion states would be induced by the hybridization 
between the orbitals. In fact, certain observed features 
can possibly be attributed to Kondo like behavior in the 
compound Ca2- x Sr x Ru0 4 (0.2 < x < 0.5)^ where the 
hybridization between orbitals is induced by the tilting 
of RuOe octahedral It is surprising that this behav- 
ior emerges from electrons which originate from bands of 
comparable width. These interesting observations natu- 
rally motivate us to explore the effect of hybridization in 
more detail. 

In this paper, we study a two-orbital Hubbard model 
with the distinct hopping integrals by combining dynam- 
ical mean field theory (DMFT) 12 i 13 i 14 i 15 with quantum 
Monte Carlo (QMC) simulations ^ML We examine the 
spin, charge and orbital fluctuations which give insight 
into the electronic properties in the regime of the OSMT. 
We further consider the effect of hybridization, which 
may be important in real materials, and show that heavy- 
fcrmion-likc behavior emerges upon introduction of the 
hybridization. The paper is organized as follows. In $U\ 
we introduce the model Hamiltonian for the two-orbital 
system and briefly explain the framework of DMFT. We 
discuss how the spin and orbital fluctuations affect the 
metal-insulator transition in flllll A brief summary is 
given in the last section. 



II. MODEL AND METHOD 

We consider the following two-orbital Hubbard Hamil- 
tonian, 
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where cj (cj aer ) creates (annihilates) an electron with 
spin o"(=|, J.) and orbital index a(= 1,2) at the ith site 
and riiaa- — c] aa .Ci aa . U (U 1 ) represents the intraband 
(interband) Coulomb interaction and J the Hund cou- 
pling. For electron hopping, we introduce 
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with the orbital-dependent nearest-neighbor hopping tff 
and the hybridization V between two orbitals. By this 
generalized model, we can study several different mod- 
els in the same framework. For V = 0, the system is 
reduced to the multi-orbital Hubbard model with the 
same (tj^ = iy) or distinct orbitals^ On the other 

(2) 

hand, for t\j — 0, the system is reduced to a correlated 
electron system coupled to localized electrons, such as 
the periodic Anderson model (J = 0) for heavy-fermion 
systemsi&iSiSSiSiiS^ or the double exchange model (J > 
0) for some transition metal oxides iSiSi 2 ^ 2 ^ For general 
choices of the parameters, we expect a variety of char- 
acteristic properties inherent in these limiting models to 
appear naturally. 

To investigate the above degenerate Hubbard model, 
we make use of DMFT ) 12 i 13 i 14 i 15 which has successfully 
been applied to various electron systems such as the sin- 
gle band Hubbard model, 27 ! 28 ! 29 ! 30 ! 31 ! 32 ! 33 ! 34 the multi- 
orbital Hubbard model, 3 ! 4 ! 6 ! 17 ! 35 ! 36 ! 37 ! 38 ! 39 ! 40 the periodic 
Anderson modeli 41 ! 42 i 43 i 44 i 45 ! 46 In DMFT, the lattice 
model is mapped to an effective impurity model, where 
local electron correlations are taken into account pre- 
cisely. The lattice Green function is then obtained via 
self-consistent conditions imposed on the impurity prob- 
lem. 

In DMFT for the multi-orbital model, the Green func- 
tion in the lattice system is given as, 



G(fc,«) =Go(M)~ -S(z). 



with 



G (k,zY 



Z + fl- £i(fc) 

-V 



-V 
z + p - e 2 (k) 



and 



VM _ ( Sn(z) S i2 (z) 



(3) 



(4) 



(5) 



where p is the chemical potential, and e a (k) is the bare 
dispersion relation for the a-th orbital. In terms of the 



density of states p(x) rescaled by the band width D a , the 
local Green function is expressed as, 
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where 



£i(z,x) = z + (i — En — D\x, 
£ 2 {z, x) = z + p - E 2 2 - D 2 x, 
v(z) = V + -E 12 (z). 
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In the following, we use the semicircular density of states 
p[x) = 2/irVl-x 2 . 

There are various numerical methods to solve the ef- 
fective impurity problem. Note that the explicit model 
Hamiltonian for the impurity system is not obtained 
straightforwardly in our case, since the lattice Green 
function has a frequency-dependent term in the non- 
diagonal element when the system has the hybridization 
V and the finite band width in both orbitals. There- 
fore, it is not necessarily most efficient to apply the ex- 
act diagonalizationSi or the two-site DMFT 28 methods 
as impurity solvers, because these methods require the 
knowledge of the explicit form of the Hamiltonian. Fur- 
thermore, self-consistent perturbation theories such as 
the iterative perturbation method and the non-crossing 
approximation are not appropriate to discuss orbital fluc- 
tuations in the vicinity of the critical point. In the 
present study, we make use of QMC to treat the impu- 
rity model at finite temperatures^ In this connection, 
we note here that the Hund coupling plays a key role in 
controlling the nature of the Mott transition in the multi- 
orbital system^ Therefore, it is important to carefully 
analyze the effect of the Hund coupling in the framework 
of QMC. To this end, we use the algorithm proposed by 
Sakai et al.fi 7 - where the Hund coupling is represented 
in terms of discrete auxiliary fields. When we solve the 
effective impurity model by means of QMC method, we 
use the Trotter time slices At = (TL)~ 4 < 1/6, where T 
is the temperature and L is the Trotter number. 

In the following, we fix the band widths as (D\, D 2 ) = 
(1.0, 2.0) and the chemical potential as p = —U/2 — U' + 
J/4 to discuss the metal-insulator transitions at half- 
filling. 



III. RESULTS 

A. Non-hybridizing bands 

Before presenting the results computed at finite tem- 
peratures, we briefly summarize the nature of the zero- 
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temperature phase diagram for V = obtained by 
DMFT together with the exact diagonalization, 3 which is 
shown in Fig. Q There are three distinct phases depend- 



where T is the time-ordering operator, n{r) — 
and r is an imaginary time. 
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FIG. 1: Phase diagram for the two-orbital Hubbard model 
with Di = 1 and D2 = 2. Note that the condition of rota- 
tional symmetry, U — U' + 2 J, is imposed (only the region of 
U > U' is relevant). In the phase (i) (phase (ii)), both bands 
are metallic (insulating), whereas in the phase (iii) the metal- 
lic state coexists with the Mott insulating state. Two lines 
along U = U' with J = and U'/U = 3/4 with J/U = 1/8 
are shown, for which thermodynamic properties at finite tem- 
peratures are examined in the text. 

ing on the strength of the interactions. It is seen that the 
metallic phase (i) remains stable up to large Coulomb 
interaction U along the line U ~ U' (small J), where 
the Mott transitions merge to a single transition. Away 
from the symmetric limit, i.e. U > U' with 2J = U — U', 
we find two separate Mott transitions in general. In be- 
tween the intermediate metallic phase (iii) appears with 
one band localized and the other itinerant. 

We now analyze the temperature dependence of the 
charge, spin and orbital fluctuations by combining 
DMFT with QMC simulations. We still restrict here to 
the case of non- hybridized bands (V = 0). Two typical 
sets of the parameters are considered, which satisfy the 
conditions (U'/U, J/U) = (3/4,1/8) and (1,0). As seen 
from Fig. ^ the Mott transitions occur at two differ- 
ent critical points U c \ ~ 3 and U C 2 ~ 4 in the former 
case, while in the latter case they are merged to a sin- 
gle Mott transition at the critical point U c ~ 7 for zero 
temperature. The charge (c), spin (s) and orbital (o) 
susceptibilities are defined as 

rP 

X~i = I drxjir), (8) 
Jo 

with 7 = c, s, o, and 

Xc (t-t') = <T|[n(r)-2][n(T')-2]>, 
Xs{t-t') = (T\ [n T (r) - n x (r)] [n T (r') -nx(r')]), 
Xo(t-t') = (T|K(r)-n 2 (r)]K(r')-n 2 (r')]),(9) 



FIG. 2: Orbital susceptibility as a function of the tempera- 
ture T for V = 0. Open (solid) symbols represent the results 
in the case U = U' and J = {& /U = 3/4 and J/U = 1/8) 
and dashed lines those for the non-interacting case. 

We first turn to the orbital fluctuations. The 
temperature-dependent orbital susceptibility is shown in 
Fig. In the non-interacting system, the orbital sus- 
ceptibility increases with decreasing temperature, and 
reaches a constant value at zero temperature. If we now 
turn on the interactions (fixing the ratios U'/U = 3/4 
and J/U = 1/8), the orbital susceptibility is suppressed 
at low temperatures. This implies that electrons tend 
to localized in each band independently such that onsite 
fluctuations are unfavorable. Eventually for U > U c \ ~ 
3, one of the orbitals is entirely localized, so that orbital 
fluctuations are suppressed completely, giving Xo — at 
T = 0. 

On the other hand, very different behavior can be seen 
along the line U 1 = U in Fig. ^ In this case, the 
orbital susceptibility is increased with growing interac- 
tions even at low temperatures. Interpreting this result 
in the context of the phase diagram in Fig. ^ we can 
say that the enhanced orbital fluctuations are relevant 
for stabilizing the metallic phase in the strong correla- 
tion regime. While such behavior is naturally expected 
for models with two equivalent orbitals, it appears even 
in systems with nonequivalent bands^ 

To examine whether the system shows metallic or in- 
sulating properties at finite temperatures, we calculate 
the charge susceptibility (compressibility). The obtained 
results are shown in Fig. [31 In the case U'/U = 3/4 
and J/U = 1/8, the system with U = 3 is located near 
the critical point between the metallic phase (i) and the 
intermediate phase (iii). With decreasing temperature 
the charge susceptibility decreases below T ~ 1. The 
appearance of a pseudogap feature in an intermediate 
temperature range gives rise to a depletion of density 
of states at the Fermi energy for both bands. Upon fur- 
ther lowering of the temperature the charge susceptibility 
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FIG. 3: Charge susceptibility as a function of the tempera- 
ture T for V = 0. Dashed lines represent the results for the 
non-interacting case. 



converges to a finite value, since the system still remains 
in a metallic phase, at least for one of the two orbitals. 
For U = 4, which corresponds to the boundary between 
the phases (ii) and (iii), the charge susceptibility at low 
temperatures is almost zero, suggesting that the system 
has become completely insulating corresponding to phase 
(ii). In contrast for U' — U we observe in an intermedi- 
ate range of U that with lowering temperature a decrease 
of the charge susceptibility is followed by an eventual in- 
crease at lowest temperatures (Fig. HJ. Comparing this 
with Fig. we see that the enhanced orbital fluctuations 
indeed have a tendency to stabilize the metallic state. 
We now move to the spin susceptibility. In Fig. 0] 




FIG. 4: The effective Curie constant XsT as a function of 
the temperatures T for U'/U = 3/4 and J/U = 1/8. Inset 
shows the results in the case U' — U and J = 0. Dashed lines 
represent the results for the non-interacting case. 

we plot the effective Curie constant XsT as a function of 
the temperature. We first look at the case of U'/U = 
3/4 and J/U = 1/8. At high temperatures, all the spin 



configurations are equally populated, so that the effective 
Curie constant takes the value 1/2 for each orbital in 
our units, yielding XsT ~ 1. When electron correlations 
are weak (U = 1), the system is still in the metallic 
phase, so that the Pauli paramagnetic behavior with a 
constant Xs emerges, leading to XsT — » as T — + 0. It 
is seen that the increase of the interactions enhances the 
spin susceptibility at low temperatures, as a result of the 
progressive trend to localize the electrons. The effective 
Curie constant is XsT = 2 when a free spin is realized in 
each orbital. It is seen that the Curie constant increases 
beyond the value of 2 with the increase of the interactions 
(U = 3,4). This means that ferromagnetic correlations 
due to the Hund coupling appear here. 

When U' = U (inset of Fig. 0}, both spin and orbital 
fluctuations are enhanced in the presence of the inter- 
actions. Accordingly, both spin and orbital susceptibili- 
ties increase at low temperatures, forming heavy-fermion 
states as far as the system stays in the metallic phase 
(see also Fig. 0). Note that for U — 6, at which the 
system is close to the Mott transition point, the spin sus- 
ceptibility is enhanced with the effective Curie constant 
X S T ~ 4/3 down to very low temperatures, as seen in 
the inset of Fig. 0] The value of 4/3 immediately follows 
if one takes into account two additional configurations 
of doubly-occupied orbital besides four magnetic config- 
urations, which are all degenerate at the metal-insulator 
transition point. Although not clearly observed in the 
temperature range shown, XsT should vanish at zero tem- 
perature for U — U' = 6, since the system is still in the 
metallic phase, as seen from Fig. ^ 

To see the above characteristic properties more clearly, 
we show the density of states for each orbital in 
Fig. |S1 which is computed by the maximum entropy 
method4£^2*22i When the interactions increase along the 
line U'/U = 3/4 and J/U = 1/8, the OSMT should oc- 
cur. Such tendency indeed appears at low temperatures 
in Fig. Ela). Although both orbitals stay in metallic 
states down to low temperatures (T = 1/6) for U = 1, 
the OSMT seems to occur for U — 2; namely one of the 
bands develops the Mott Hubbard gap, while the other 
band still remains metallic. At a first glance, this result 
is slightly different from the ground-state phase diagram 
shown in Fig. ^ where the system is in the phase (i) 
even at U — 2. However, this deviation is naturally un- 
derstood if we take into account the fact that for U = 2, 
the narrower band is already in a highly correlated metal- 
lic state, so that the sharp quasi-particle peak immedi- 
ately disappears as the temperature increases beyond the 
small characteristic energy scale. This explains the be- 
havior observed in the density of states at T = 1/6. For 
U = 3, both bands are insulating at T = 1/6 (the system 
is almost on the boundary between the phases (ii) and 
(iii) at zero temperature). 

In the case U' — U, as expected we encounter the qual- 
itatively different behavior shown in Fig. [SJ In this case, 
both bands gradually develop quasi-particle peaks as the 
interactions increase, and they still remain metallic even 
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FIG. 5: Density of states for the degenerate Hubbard model 
(Di,Z?2) = (1.0, 2.0). The data are for the temperatures T = 
2, 1, 1/2 and 1/6 from the top to the bottom. 




FIG. 6: Solid (dashed) lines represent the density of states for 
the orbital a = 1 (a = 2) when (Di,D 2 ) = (1.0, 2.0) at T = 
1/6 with the fixed parameters of U'/U = 3/4 and J/U = 1/8. 
The data are plotted for V = 0.0, 0.25, 0.5, 0.75, 1.0, 1.25 and 
1.5 from top to bottom. 



at U = U' = 3. As mentioned above, all these features 
which arc in contrast to the situation for U' ^ U, are 
caused by the special symmetry for U — U' , which gives 
rise to equally enhanced spin and orbital fluctuations. 

B. Hybridization between distinct orbitals 

We have so far treated the degenerate Hubbard model, 
in which two types of orbitals do not mix with each other. 
In our treatment with DMFT, the Mott insulating phase 
(ii) as well as the intermediate phase (iii) may be unstable 
against certain perturbations. There may be several pos- 
sible mechanisms that stabilize such insulating phases. 
One of the mechanisms, which may play an important 
role in real materials, is the hybridization between the 
two distinct orbitals. We address the effect in this sec- 
tion. 

This hybridization effect is relevant in some real ma- 
terials. For instance, in the compound Ca2- x Sr x RuC>4, 7 
the hybridization between {a, /3} and 7 orbitals is in- 
duced by the tilting of RuOe octahedra in the region 
of Ca-doping 0.2 < x < 0.5, 11 . This leads to Kondo- 
lattice like effective model and may be connected with 



the reported heavy fermion behavior^ similar to some 
/-electron systems. This interesting aspect motivates 
us to study the mixing effect between the localized and 
itinerant electrons in the intermediate phase (iii). More- 
over the compound La n +iNi n 03 n +i— possesses hybridiza- 
tion between d 3z 2_ r 2 and d x i_ y i orbitals in the e g sub- 
shell. The OSMT may lead to the metallic but the less- 
conducting state is realized below the critical tempera- 
ture T c = Consequently we would like also to 
explore how the hybridization of different-type <i-bands 
affects electronic properties especially around the OSMT. 

We study the general case with U' 7^ U and J 7^ in 
the presence of the hybridization V. In Fig. |B| the den- 
sity of states calculated by the maximum entropy method 
is shown for different choices of V. We start with the 
weak coupling case, {7 = 1, where the metallic states are 
realized in both orbitals at V = 0. Although the intro- 
duction of small V does not alter the nature of the ground 
state, further increase of V splits the density of states 
(V = 1.5), signaling the formation of the band insulator: 
namely all kinds of elementary excitations possess the 
gap. In contrast, we encounter different behavior when 
electron interactions are increased up to U = 2 and 3. 
In these parameters, the system at V = shows the in- 
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FIG. 7: Charge, spin and orbital susceptibilities as a function 
of the hybridization V at the temperature T = 1/6. 



termediate or Mott-insulating properties at T = 1/6. It 
is seen that the density of states around the Fermi level 
increases as V increases. For U = 2, the intermediate 
state is first changed to a metallic state, where the quasi- 
particle peaks appear in both orbitals (V = 0.75,1.0). 
For fairly large V, both bands fall into the renormalized 
band insulator (V = 1.5). Similarly, for [7 = 3, the hy- 
bridization first drives the Mott-insulating state to an 
intermediate one, as is clearly seen at V = 0.75, which is 
followed by two successive transitions as is the case for 
U = 2. 

The above characteristic properties also emerge in the 
charge, spin and orbital susceptibilities at low tempera- 
ture, as shown in Fig. [7| For weak interactions (U = 1), 
the charge susceptibility \c monotonically decreases with 
the increase of V. When electron correlations become 
strong, the non-monotonic behavior appears in Xc- the 
charge fluctuations, which are suppressed at V — 0, are 
somewhat recovered by the hybridization, which leads to 
metallic behavior. For large V, Xc is again suppressed 
since the system turns into a band insulator. We can see 
that the orbital susceptibility exhibits non-monotonic be- 
havior similar to the charge susceptibility, the origin of 
which is essentially the same as in Xc\ the orbital fluctu- 



FIG. 8: (a) Effective Curie constant as a function of the 
temperature and (b) density of states in the narrower band 
(a = 1) at T — 1/4 for an extreme choice of the band- 
widths, (Di,D 2 ) = (1.0, 10.0). The density of states for the 
wider band is not shown here. The other parameters are 
U = 4.0, U' = 3.0 and J = 0.5. 



ations suppressed at V = are recovered by V, and then 
the formation of the band insulator causes the gradual 
decrease of Xo- In contrast, the spin susceptibility mono- 
tonically decreases with the increase of V irrespective of 
the strength of the interactions. As discussed for V = 0, 
the effective spin is enhanced by ferromagnetic fluctu- 
ations due to the Hund coupling in the insulating and 
intermediate phases. Upon introducing the hybridiza- 
tion in these phases, the ferromagnetic fluctuations are 
suppressed, leading to the monotonic decrease of the ef- 
fective Curie constant. 

From the above observations, we can say that the in- 
troduction of appropriate hybridization induces heavy- 
fermion metallic behavior. In fact, this tendency can 
be observed more clearly in an extreme choice of the 
bandwidths, (D 1 ,D 2 ) = (1.0,10.0), shown in Fig. 
At V = 0.0, the system is in the intermediate phase, so 
that the completely localized states [Fig. [S] (b)] appear in 
the narrower band in the background of the nearly free 
bands. This double structure in the system gives rise to 
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FIG. 9: Schematic phase diagram for the two-orbital Hub- 
bard model with finite hybridization between two orbitals. 
Solid lines represent the phase boundaries between the metal- 
lic and insulating phases. Dashed line indicates the crossover 
between the Mott insulator and the Kondo insulator. 



two peaks in the temperature-dependent effective Curie 
constant, as shown in Fig. |8](a). Since the completely 
localized state plays a role of the /-state in the Anderson 
lattice models a " heavy- fermion" peak appears at the 
Fermi energy in the presence of V , which is essentially 
the same as that observed in Fig. |SJ 

Finally, some comments are in order on the phase di- 
agram at zero temperature. In our approach, it is not 
easy to deal with the system at very low temperatures, 
since QMC simulations suffer from minus sign problems. 
Nevertheless, we may give some qualitative arguments 
on the expected phase diagram at zero temperature. As 
discussed above, the metallic phase (i) is not so sensi- 
tive to V as far as it is small. This is also the case 
for the completely insulating phase (ii). In contrast, a 
more subtle situation appears in the intermediate phase 
(iii). As mentioned above, the intermediate phase ex- 
hibits Kondo-like heavy fermion behavior at low temper- 
atures in the presence of V. Recall, however, that we 
are now concerned with the half-filled band. Therefore, 
this Kondo-like metallic phase should acquire a Kondo- 
insulating gap due to commensurability at zero tempera- 
ture. We would thus say that the intermediate phase (iii) 
is changed into the Kondo-insulator with a tiny excitation 
gap in the presence of V at zero temperature. Accord- 
ingly, the sharp transition between the phases (ii) and 
(iii) at V — may be smeared and changed to crossover 
behavior. These considerations lead us to a schematic de- 
scription of phase diagram for the two-orbital model with 
mixing between the distinct orbitals, as shown in Fig. 
On the line of V — 0, the OSMT, which may occur in 
general choices of the parameters, separates the phase at 
V = into three regions. The metallic phase for small 
U is simply driven to the band-insulator (iv) beyond a 
certain critical value of hybridization. The intermediate 



phase at V = is changed to the Kondo-insulator in the 
presence of any finite V. This insulating state first under- 
goes a phase transition to the metallic phase, and eventu- 
ally enters the band-insulator as V increases. The com- 
pletely Mott insulating phase first shows a crossover to 
the Kondo insulator, which is further driven to the metal- 
lic phase and then to the band-insulating phase. Note 
that at finite temperatures above the Kondo-insulating 
gap, we can observe a Kondo- type heavy fermion behav- 
ior in the intermediate phase with finite V. 



IV. SUMMARY 

We have investigated the degenerate Hubbard model 
with distinct hopping integrals by combining DMFT with 
QMC simulations. By examining the spin, charge and or- 
bital susceptibilities calculated at finite temperatures, we 
have clarified that equally enhanced spin and orbital fluc- 
tuations play a vital role on stabilizing the metallic states 
in the multi-orbital systems. This remarkable effect is 
responsible for whether the system undergoes a single 
Mott transition or OSMTs. Also, we have discussed how 
the phase diagram at finite temperatures slightly deviates 
from the ground-state one because of smearing effect of 
the narrow quasi-particle peak. 

We have further explored the effect of the hybridiza- 
tion between the distinct orbitals, and have found that it 
plays a crucial role especially around the OSMT. The in- 
troduction of the hybridization in the intermediate phase 
enhances the charge and orbital fluctuations, inducing 
the metallic phase with a sharp quasi-particle peak. Ac- 
cordingly, Kondo-like heavy fermion states show up at 
finite temperatures, which eventually drop in the Kondo 
insulating phase for our half-filled bands. We have also 
pointed out that the hybridization effect smears the sharp 
OSMT at zero temperature, and changes it to a crossover 
behavior. Nevertheless, we can still observe the OSMT 
at finite temperatures. 

In this paper, we have used QMC as an impurity solver 
in DMFT, which is not powerful enough to treat proper- 
ties at very low temperatures. Therefore, it is desirable 
to exploit a complementary approach to study such low- 
temperature properties more precisely, although we have 
arrived at a reasonable phase diagram at zero tempera- 
ture. Various remaining open problems could not be ad- 
dressed in the present study. One of the most important 
issues to explore is magnetism of the system, which has 
not been seen here, since we have restricted our atten- 
tion to the paramagnetic phase. This problem is under 
consideration. 
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